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In the previous paper, we studied the random-mass Dirac fermion in one dimension by using 
the transfer-matrix methods. We furthermore employed the imaginary vector potential methods 
for calculating the localization lengths. Especially we investigated effects of the nonlocal but short- 
range correlations of the random mass. In this paper, we shall study effects of the long-range 
correlations of the random mass especially on the delocalization transition and singular behaviours 
at the band center. We calculate localization lengths and density of states for various nonlocally 
correlated random mass. We show that there occurs a "phase transition" as the correlation length 
of the random Dirac mass is varied. The Thouless formula, which relates the density of states and 
the localization lengths, plays an important role in our investigation. 



I. INTRODUCTION 

In the previous papers j^J|], we studied the random-mass Dirac fermions in one dimension by using the transfer- 
matrix methods (TMM) and imaginary- vector potential methods (IVPM). We calculated the density of states, the 
typical and mean localization lengths, and the multifractal scalings as a function of the energy and correlation length 
of the random mass. The results are in good agreement with the available analytical calculations ||,^. We also 
obtained the relation between the correlation length of the random mass and the typical locahzation length for the 
short-range correlations ||^. 

In this paper, we shall study effects of the long-range correlations, especially on the delocalization transition and 
singular behaviours at the band center. This problem is interesting for various reasons. Very recently finite mobility 
edges are observed in the one-dimensional Anderson model with the long-range correlated potentials though it is 
widely believed that (almost) all states are localized in spatial one dimension [^j6|. Similar phenomena are observed 
also in the aperiodic Kronig-Penney model 0. In contrast to the above models, the system of the random-mass 
Dirac fermions belongs to the universality class of the chiral orthogonal ensemble and therefore the extended states 
exist at the band center though the others are localized for the white-noise random mass. Long-range correlations of 
the disorder may change this feature. This problem has not been addressed in detail so far. As the model is closely 
related with the random-bond XY model, the random Ising model, etc., the results are interesting both theoretically 
and experimentally. 

This paper is organized as follows. The model and numerical methods are explained in the previous paper |2|. In 
Sec. 2, we shall explain how to make a long-range correlated random Dirac mass numerically. We shall consider two 
types of the telegraphic random mass, i.e., in the first one magnitudes of the random mass are long-range correlated 
random variables (LRCRV) with fixed interval distances between kinks whereas in the second one interval distances 
between kinks are LRCRV with a fixed random mass magnitude. In the present study, we employ the first one in 
the above for it generates various correlations of the Dirac mass more easily than the second one. In Sec. 3, numerical 
results are given. We calculate the localization lengths and the density of states (DOS) directly by using the TMM 
and IVPM though in most of studies of one-dimensional random models elaborated techniques like the Hamiltonian 
mapping, the renormalization group, etc. are used for calculating the Lyapunov exponent. For the white-noise and 
short-range correlated random mass, the DOS p{E) diverges as p{E) cx ^| jogg|3 where E is the energy measured 
from the band center. Similarly the typical localization length £,{E) cx | \ogE\. This means that there are extended 
states at the band center. We study how these behaviours change by the existence of the long-range correlations of 
the Dirac mass. We show that there is a kind of "phase transition" as the correlation length is varied. Section 4 is 
devoted to conclusion. 



II. MODELS AND LONG-RANGE CORRELATED RANDOM MASS 



Hamiltonian of the random-mass Dirac fermion is given by 
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h = -ia'd.^+m{x)<7y, (1) 

where a are the PauU matrices and m(x) is the telegraphic random mass. It is known that almost all energy eigenstates 
of the Hamiltonian (|l|) away from the band center are localized for short-range correlated m{x) |^ . There we calculated 
localization length of each state by introducing an imaginary- vector potential (IVP) g as suggested by the work of 
Hatano and Nelson . In an IVP, all localized states tend to extend and a "critical" value of IVP gc at which a state 
becomes extended, "determines" its localization length a.t g — 

In this paper we shall consider the long-range correlated configurations of m(x) and study delocalized states by the 
TMM and IVPM. Numerical methods which generate potentials correlated by the power law was recently invented 
by Izrailev et.al. and Herbut We briefly review it for we shall use it for the present studies. 

Let us consider a spatial lattice and random potential e{n) sitting on the sites which has correlation x(n), 

{ e{m) e{n) ]cns ^ C x{\m - n\), (2) 

where n and m are site indices, and C is a constant. Here we suppose that xl"-) is given and define c(fc) from xl"-) 
as follows, 

c{k) = {my\ (3) 

where x(^) the Fourier transformation of x(^); 



oo 



X(fc)= J2 xWe''^". (4) 

n— — oo 

We also introduce another random site potential r{n) which is distributed imiformly in the range [—1,1] with the 
delta-function white-noise correlation; 

[ r{m) r{n) ]ons oc Smn- (5) 

Then we can construct random potential e{n) as follows by using the white-noise potential r(m) and the inverse 
Fourier transformation of c(fc), 

oo 

e(n) = C r(m) c(m — n), (6) 

m— — oo 

where C is another constant. However, we cannot perform the infinite summation in in the practical numerical 
calculation, and therefore we replace it with the finite summation as follows, 

k/2 

e{n) — C r{m)c{m — n). (7) 

m=-fc/2 

where k is a finite number which is much larger than the number of sites in the system. We use the periodic boundary 
condition for the numerical calculation in subsequent sections and set k ~ 10"* -I- (number of sites). 
We can easily verify that the random variable e(n) in satisfies Eq.(|^). If we choose c(m) as 

with some constant a, then the correlator of the random potential becomes 

[ ]ons ~ 1 ^|2a-l - (9) 

With this procedure, we can generate power-law correlated random potential. Similarly, we can generate exponentially 
correlated random potential if we choose the modified Bessel function Kq as c(m). 
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In the practical calculation, there is a subtle point in the normalization of the random potential e(n) in Eq.(||) when 
we vary the system size L |[ic| ]. We determine the normalization of e(m) (namely the value of C ) in such a way that 
the typical value of e(m) does not change. Then as a result, the strength of correlation does not change for various 
system size L if we fix |m — n| in Eq.(^. 

The above method generates the random potential at each site. Therefore, when we fix the distances between kinks 
and vary the magnitudes of m{x) for generating random m{x)^ we can directly use the above methods (see Fig.l). We 
can easily extend TMM to the above type of random mass |TO(a:)|. Numerical studies of this system will be reported 
in the following section. 

m(x) 




FIG. 1. An example of m{x) 

On the other hand, we can also use the above method for generating random distances between kinks with fixed 
value of |TO(a;)|. This is the subject in our previous papers 0,^. 



III. LONG-RANGE CORRELATED DISORDERS 



We focus on the system with exponential and power-law correlated random mass in this section. Each correlation 
is parameterized as follows; 



m{x) m{y) 
m{x) m(y) 
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FIG. 2. An example of m{x) in the case of white-noise disorder(top left), exponentially correlated (top right) and power-law 
correlation(bottom): We set 2000 kinks in these systems. By the numerical analysis, correlation length (per system size) A is 
calculated as 0.16 in exponentially correlated case and is 0.68 in power-law case. 
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where g,X,C and Upw are constants. Dimensions of the parameters are [A] = M~^, [g] = M^, [C] = M^"""™, 
respectively. As we explained in Sec. 2, we consider the random mass of random magnitude with fixed distances 
between kinks. Prototype of the configuration is given in Fig. 2. In this case m{x) is generated rather similarly to 
e{n) in Sec. 2 though the continuum space is considered here instead of the spatial lattice. As in the previous paper 
1^, we use the IVPM and TMM in order to calculate the localization lengths. 

We first calculate the correlation [m(x)m(0)]cns numerically which is expected to exhibit the exponential and 
power-law decay, respectively. The result is shown in Fig. 3. The correlator actually shows the expected behaviour 
in each case, though the parameters of the correlation are slightly smaller than the analytical values, which are A in 
exponential-decay case and ap^ in power-law decay case, respectively. We think that this is due to the finiteness of 
numbers of kinks and/or finite Fourier transformation in the scheme which we explained in the previous section. (The 
values of apw , A and the parameters in Tables 1 and 2 are calculated directly from the two-point correlation function 
of the random mass obtained numerically.) 
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FIG. 3. Correlation of m{x) in the case of exponentially correlation (left) and power-law correlation (right): Each case cor- 
responds to random m{x) shown in Fig. 2. In these cases, we can consider that randomness has exponential or power-law 
correlation for about 10 percent of system size. We fix the parameters in Table 1 (exponential case) and 2 (power-law case) 
using this numerical result. (In order to fix these parameters, we used the data points in the range 0.003 < a < 0.015 for 
exponential case (0.01 < a < 0.05 for gX — 12) and 0.01 < a < 0.1 for power-law case.) 



Next we show the energy dependence of typical localization length in the system with the white-noise random mass. 
(We obtain "typical" localization length as a function of energy E by averaging localization lengths of the eigenstates 
within a small range of energy AE.) In Fig. 4, we show the ratio of the numerical results to the available analytical 
expression which takes the form 

C(i?)(x|log(i?/2g)|, (11) 

where the parameter g gives the magnitude of m{x), 

[ m{x) m{y) ]cns ^ g S{x ~ y). (12) 

It is easily seen that the limit A ^ in Eq.(|o|) corresponds to Eq.(p^). As Fig. 4 shows, our numerical calculations 
and the analytical expression ( [Tl| ) are in good agreement. This guarantees validity of our methods of calculation, the 
TMM and IVPM. 
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FIG. 4. The energy dependence of localization length in the case of white-noise disorder : We set L(system size)=50, 1000 
kinks in the system and 5 = 1. Ratio is normalized at E — 0.10 and data are averaged within small energy slice AE, We set 
AE = 0.02 here. 
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We calculated the energy dependences of the localization length in the case of exponentially correlated disorder in 
Eq.(|lO|). In Fig. 5, we show the ratio of numerical calculations to the analytical expression of the white-noise case 
I \og{E/2g)\. We vary the correlation length A of m{x). The ratio is almost constant for small gX. However as gX is 
getting large, the numerical results obviously deviate from | \og{E/2g)\. 
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FIG. 5. The energy dependence of localization length in the case of exponentially correlated random mass: We set L(system 
size)= 50 and 1000 kinks in the system. The ratios are normalized a,t E — 0.1 and data are averaged within energy slice 
AE = 0.02. Values of A and g are given in Table 1. The ratio is almost constant for small gX, but it deviates from constant in 
the case of large gX. 



Table 1. Parameters of exponential correlation. 
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We show the energy dependence of the localization length in the case gX = 12 in Figs. 5 (bottom) and 6. We can 
conclude that the localization length diverges only a.t E — even for very large gX. This means that exponential 
correlation of the random mass gives no significant effect on the Anderson transition. Actually this result can be 
expected from the study in the previous paper pi. 




FIG. 6. The energy dependence of localization length in the case of exponentially correlated random mass (especially the 
case of large g\): This data corresponds to the one at the bottom in Fig. 5. Data are averaged within energy slice AE = 0.02. 
Divergence of localization length is found only at the point E = 0. 
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However the result also shows that singular behaviour of the localization length at E = for the exponentially 
correlated m{x) (with large gX) may be different from that of the white-noise case. From the analytical study by using 
supersymmetry the localization length is obtained in powers of gX and more dominant terms like \ \og{E /2g)\^ 
etc. appear in higher-order terms of gX. 

We therefore investigate the behaviour of the localization length near the band center E — Q rather in detail. In 
Fig. 7 we show the ratio of the numerically obtained localization length to | \og{E /2g)\^ ^ where 5 takes the value 1,2,3 
and 4. In some energy regions, the energy dependence is fitted better by the power 5 — 2 ov Z than 5=1. 




Energy Energy 



FIG. 7. Ratio of the Localization length to | \og[E /2g)\^ : Here we set L(system size)= 50, and 5 takes the value 1,2,3 and 4. 
Data are averaged within energy slice A_B = 0.01. 

We can expect the above result by the analytical study. In the previous paper, Q the DOS p{E) is obtained 
analytically as follows for the exponentially correlated random mass; 

(E) = A2 ^3 , y,. 

' f |log(£;/2g)|3 ||log(i?/25)|4 f^\\og{E /2g)r ^ ^ 

where Ai (i = 1, 2, 3) are polynomials of gX, 

A, = l-^{gXf-—{gXf + ---, 

2 15^^ ' 105^-^ ' 
A2 = -3(.9A) - 2{gXf + • • • , 

A3 = AigXf + ■■■. (14) 

In order to verify Eq.(^3|), we calculate the DOS numerically for relatively small gX{< 1) and compare the result 
with Eq.(^3|). From Fig. 8, it is obvious that p{E) up to the third-order term of {gX) is in better agreement with the 
numerical result than the first-order one. 

Using Thouless formula for ID random hopping tight binding model (whose low-energy field theory is the 
random-mass Dirac fermion), 

1 /-l^l ,E 



oc 



\og{—)p{E')dE\ (15) 



we can show that corresponding to the second and the third term in ( p^ ) there are contributions to the localization 
length like | \og{E /2g)\'^ and | \og{E /2g)\'^ if these terms become dominant in the DOS. It is also expected that the ratio 
is not constant for any range of energy in the case 6 = 4, because DOS has no term proportional to 1/E\ log{E/2g)\^. 

These numerical results and the above discussion indicate that the nonlocal correlations of the random mass generate 
nontrivial effect on the localization. One may expect that for large gX the perturbative calculation in powers of gX 
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breaks down for the DOS and the localization length and p{E) and ^{E) exhibit different singular behaviours at 
E = 0. Especially the terms like | log-Ej"^ in ^{E) may be summed up and the localization length may behave as 



(16) 



for large gX where /3 is a constant. In Fig. 9 we show the log-log plot of (,{E) v.s. energy. The numerical data are 
obviously on a straight line and therefore the singular behaviour ([l6| ) is verified. 

From the Thouless formula (|l5|), we also expect the power behaviour of the DOS. Actually it is not so difficult to 
show that for the DOS parameterized as 

p{E)^E-^\\ogE\-'', (17) 
where rj and 7 are parameters, the Thouless formula ([l5| ) gives the following localization length, 

^{E) - E-^+'^\logE\-^+'^, for Tjj^l, 

^{E) \\ogE\~^+^, for 77 = 1. (18) 
In Fig. 10 we show the calculations of the DOS for gX = 12. It is obvious that the DOS p{E) behaves as 

p{E) ^ E-^' . (19) 

We cannot fix the exponent 7 in Eqs. ( ^jj ) and (|l^) from these numerical calculations. However, we can conclude 
that T] ^ 1 in the case gX = 12. This confirms the power-law behaviour of the localization length near the band center 
E = 0, which differs from the behaviour of the white-noise disorder case. 
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FIG. 8. The energy dependence of DOS in the case of exponentially correlated random mass (especially the case of large 
gX): We show the ratio of DOS we calculated numerically to the one obtained analytically. Ratio is normalized at i? = 0.10 
and data are averaged within energy slice AE = 0.02. We set gX = 0.80 here. (This corresponds to the case of the top left 
figure in Fig. 5.) In order to show the effect of exponential correlation, we use DOS obtained analytically up to 0th and 3rd 
order of gX expansion, respectively. 




Energy (Log) 

FIG. 9. The energy dependence of localization length in the case of exponentially correlated random mass (especially the 
case of large gX): Data are averaged within energy slice AE = 0.005. We set gX = 12 here. (This corresponds to the case of 
the bottom figure in Fig. 5.) x^{P^^ freedom) value of this fitting is 0.36. 
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FIG. 10. The energy dependence of DOS in the case of exponentially correlated random mass (especially the case of large 
gX): Data are averaged within energy slice AE = 0.005. We set gX = 12 here. (This corresponds to the case of the bottom 
figure in Fig.5.) We set energy slice AE = 0.002,0.006 and 0.01 for E < 0.005,0.005 < E < 0.01 and 0.01 < E respectively. 
The interval of data points are 5 x 10"*, 1 x 10"^ and 2.5 x 10"^ for E < 0.015, 0.015 < -B < 0.02 and 0.02 < E respectively. 
DOS is normalized at log E — —4.0. value of this fitting is 0.64. 



Let us turn to the study on the power-law correlated random mass in Eq.(|10| 

C 
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Here we emphasize significance of the power-law correlation of the randomness to the Anderson transition etc. This 
is due to the scale-free property of the power-law correlation, and this type of correlations may not be negligible even 
though system size is very large. 
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FIG. 11. The energy dependence of localization length in the case of power-law correlated random mass: We set L(system 
size)= 50 and 1000 kinks in the system, a^m's and C's in Eq.(p^ are given in Table 2. Localization length is averaged within 
energy slice AE = 0.003, 0.004, 0.006 and 0.01 for < £ < 0.01, 0, 01 < S < 0.015, 0.015 < -E < 0.02 and 0.02 < E respectively. 
The interval of data points are 5 x 10"*, 1 x 10"^ and 2.5 x 10"^ for E < 0.015, 0.015 < -B < 0.02 and 0.02 < E respectively. 
Divergences of localization length are observed only at E — 0. 
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In Fig.ll, we show the energy dependences of the localization length in the system with power-law correlated 
random mass. We choose value of the parameter apw as in Table 2. These results show that localization length 
diverges only at E' = 0, as in the case of white-noise and exponentially correlated disorder. This behaviour does not 
depend on the value of apw 

We also study the system size dependence of the localization length. Results are shown in Fig. 12. We conclude 
that extended states exist only at the band center E = 0. 
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FIG. 12. The system size dependence of typical localization length in the case of power-law correlation: Parameters except 
for the system size and number of kinks are the same used in the top right figure in Fig.ll. (Here we fix the number of kinks 
per system size.) Divergence of localization length is observed at £ = in each case. 





apw=0.68 






-0.067x-0.21 





Energy (Log) 





0.45 


N 


0.4 


'J: 




em 


0.35 


yst 


03 


in 






0.25 


J 




J 


0.2 


■3 






0.15 


'a. 






0.1 



apw=0.34 



.1 -1.2 



Energy (Log) 



apw=0.099 



3 

■5, 0.15 ^ -0.080x-0.19 

H 0.1 ■ 

0.05 ■ 







apw=0.6S 


-0.53x -4.72 






5 -6 


.5 -5 -4.5 

Energy (Log) 


-4 -3 






apw=0.34 


-0.51x-4.17 






5 -6 


.5 -5 -4.5 

Energy (Log) 


-4 -3 





apw=0.099 


-0.53x - 4.22 r+ 


S 



-6.5 



-5.5 



Energy (Log) 



Energy (Log) 



FIG. 13. Log-uni and log-log plot of the data shown in Fig. 9: The left column is log-uni and the right one is log-log plot. 
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Table 3. Results of linear fit to log-uni and log-log plot: We use the data in the range —6 < logE < —3.5 (—6 < logi5 < —4.4 
for Qpu, = 0.68). The incline in the table corresponds to —fJ in Eq.(^). values for log-log plot are estimated smaller than 
log-uni plot. 
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In order to investigate the singular behaviours of the Dirac fermions with the long-range correlated random mass, 
we change the linear scale of energy and localization length in Fig. 11 to log-uni and log-log scale plots in Fig. 13. 
If the data in the log-uni scale are on a straight line, the energy dependence of localization length £,{E) is given by 



^{E) (x\logE\+ const. (E-^O), 



(20) 



as in the white-noise case. On the other side, if the data in the log-log scale are on a straight line, the energy 
dependence of £,{E) is given as Eq.([l6|) £,{E) oc E~^ up to possible | logi?| corrections. From Fig. 13 and Table 3, we 
can conclude that -E-dependence of ^{E) is given by Eq.jl^). The exponent f3 is estimated as in Table 3. 

The above singular behaviour in the present system a.t E — 0, i.e., Eq.([T^) is also confirmed by the calculation of 
the DOS near the band center E — 0. We show the energy dependence of the DOS in Fig. 14. 
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FIG. 14. The energy dependence of DOS in the case of power-law correlated random mass: DOS is calculated as the number 
of states within energy slice AE = 0.002, 0.006 and 0.01 for E < 0.005, 0.005 < -B < 0.01 and 0.01 < E respectively. The 
interval of data points are 5 x 10"'', 1 x 10"^ and 2.5 x 10"^ for E < 0.015, 0.015 < E < 0.02 and 0.02 < E respectively. DOS 
are normalized at log E = —4.0. 
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Table 4. Results of linear fit to DOS log-log plot: We use the data in the range —6 < log-E < —3.5 (—6 < log.B < —4.4 for 
ctpw = 0.68). The incline in the table corresponds to —f3' in Eq.(pl|). 



The result indicates that the DOS p{E) diverges as 

p{E) oc E-^' {E 0), (21) 

with the constant /?' given in Table 4. From the Thoulcss formida Eq. (p^) , f] in Eq.(p^) and (3' in Eq.(pll) must satisfy 
the relation like 

f3 + f3' = l. (22) 

From Tables 3 and 4, it is seen that P + P' is almost unity as it is expected. In order to determine the exponent 7 of 
the possible | logi?| correction in Eqs.(p7|) and (p^), calculation of the localization and the DOS for wider range of E 
is required. 
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In the above we concluded that the power-law correlation of the random mass in Eq. ( |lO| ) does not influence the phase 
structure of the system and extended states appear only at E = 0. This result may seem to be in contradiction to the 
results obtained in the previous papers in which nonlocally-correlated Anderson, random hopping and random 

Kronig-Penney (random delta potential) models are studied and in each model the presence of a nontrivial mobility 
edge is observed at a finite energy. However, the correlations of randomness in their models are different from those 
of the present study. Actually in one of these papers , the following correlation was used for the one-dimensional 
random hopping tight binding(RHTB) model, 

[(e,-e,f],„, K |i-jr-\ (23) 

where w is a non-negative constant and they concluded that a delocalized phase appears in a finite region around 
i? = in the case of w = 1.5. As the random- mass Dirac fermion is a low-energy effective field theory of the RHTB 
model, the above result is relevant to the present study. As in the present study we consider the case a; < 1, their 
results are consistent with ours. 

We use the random variables in Eq.(|2^) for m(x) in the present field-theory model and calculate the localization 
length as before. The result of the numerical study is shown in Fig. 15. From these calculations, we cannot find any 
finite region of extended states around E — even in the case lo — 1.5. However, the system size may not be large 
enough for this calculation. 
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FIG. 15. Localization length and energy of the eigenstates in correlated disorder system: Correlation of m{x) is given in 
Eq.(E3). We set L(system size)=50 and 200 kinks in the system. We set uj — At the left and lo = 1.5 at the right. 



IV. CONCLUSION 



In this paper, we studied the Dirac fermion with the long-range correlated random mass by using the TMM and 
IVPM. We are especially interested in the delocalization transition or the existence of a finite mobility edge. Recent 
studies on the Anderson model, random hopping model and random Kronig-Penney model indicate the existence of 
finite mobility edges. However we found that extended states exist only at the band center i? = as in the white-noise 
random variable case. 

We also studied the localization length as a function of E rather in detail. We found that in the case of short-range 
correlation £,{E) cx |logi?| for small E whereas ^{E) oc E^^ (/3 ^ 0.5) for the long-range correlation of the random 
mass. The above conclusion is supported by the calculation of the DOS through the Thouless formula. We hope 
that the above behaviour of the localization length is observed by experiment of the random spin model because the 
localization length is directly related with the spin-spin correlation length there. 

We can also calculate the multi-fractal scaling indices directly from wave functions of the random-mass Dirac 
fermions. Values of the indices may depend on the decay power of the correlation of the random mass. Results will 
be reported in a future publication |l3| . 
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